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The Josephson effects in two weakly linked Bose- 
Einstein condensates have been studied recently. In 
this letter, we study the equations derived by Gio- 
vanazzi et. al. [Phys. Rev. Lett. 84, 4521 (2000)] 
focusing on the effects of the initial acceleration 
and the velocity of the barrier on the "dc" current. 
We find that the dc current has lifetime which criti- 
cally depends on the moving velocity of the barrier. 
Moreover, the influence of the initial acceleration is 
also investigated and found to be crucial for the ex- 
perimental realization of the effects. 
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The Josephson effects ( JE's) as a paradigm of 
the phase coherence manifestation in a macro- 
scopic quantum system, have been observed 
in superconductors [jjj, and demonstrated in 
two weakly linked supcrfiuid 3 He-B reservoirs 
g. Since magnetic and optical traps can be 
tailored and biased with high accuracy || ||, 
the weakly interacting Bose-Einstein conden- 
sate (BEC) can provide a further context for 
JE's and reveal novel properties that might 
not be accessible with other systems. Re- 
cently the dc and ac Josephson effects in two 
weakly linked BECs have been extensively stud- 
ied |]]^]. These authors suggested that as the 
barrier between the two trapped BECs moves 
adiabatically across the trapping potential, a dc 
current of atoms between two condensates can 
be found. In analog of the voltage-current char- 
acteristic in superconducting Josephson junc- 
tion (SJJ), there exists a critical velocity of the 
barrier, at which an abrupt transition from the 
dc to ac current occurs. In this letter, we study 
the model introduced in Ref. || focusing on the 
effects of the initial acceleration and the veloc- 



ity of the barrier on the "dc" current. We find 
that the dc current has lifetime which depends 
on the velocity of the barrier and is sensitive 
on the choice of the initial conditions. To con- 
sider the experimental observability of this phe- 
nomenon, we investigate the influence of the 
initial acceleration and find it plays a crucial 
role. 

The interacting BECs in a trap at zero tem- 
perature can be described by a macroscopic 
wave function ^(r,i), having the meaning of 
an order parameter and satisfying the Gross- 
Pitaevskii equation 



where Hq = 



[ff +.g|*(r,i)| 2 Mr,i), (1) 
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V ext (r,t) and g 



AwTi a/m with m the atomic mass, and a the 
s-wave scattering length of the atoms. Consid- 
ering the system proposed in Ref. || , a double- 
well trap produced by a far off-resonance laser 
barrier, Vi aser (z) = V Q exp[-(z - Vf/X 2 ), which 
cuts a single trapped condensate into two parts 
0. So, the external potential is given by the 
magnetic trap and the laser barrier V ex t (r, t) — 

Vtrap{r) + Vl 

a scr 

CM)- 

By solving variationally the GPE using the 
ansatz: \- (r,i) = -i/>i(i)0i(r-)+?/'2(O ( / , 2(r), where 
■01.2 = \J N\,2{t)e %ei ' 2 ^ are complex time- 
dependent amplitudes, Ni^(t) and 6^2 are the 
number of atoms and the phase of the two con- 
densates respectively. The trial wave functions 
cf>i 2 ( r ) are orthonormal and can be interpreted 
as approximate ground state solutions of the 
GPE of two well respectively. Then the equa- 
tions of the motion for the relative population 
p(t) = (N 1 (t)-N 2 (t))/N and phase 9 = 9 2 - 61 
between the two condensates should be ||U 



= -Fl 



Kp 



r cos 9 + Cp, 



(2) 
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where K = - § / dr0 x [# + gN(j>\\<h, C = 
ag-J^dr and F = - 0|]f Vi aser 

. These equations describe the dynamics of the 
Bose Josephson junction (BJJ). For the con- 
venience we set the parameters as the same 
as in Ref. §, K = 4.82 x lO^ms" 1 , C = 
1.23 ms -1 and F = 1.06 ms _1 /Ltm _1 . 

The Josephson current is defined as the cur- 
rent of atoms across the barrier, which is J — 
N% = —Ni- Considering the total number of 
atoms is constant , we introduce the normal- 
ized current j = J/N . From the Eq. (||), we 
have 



j= 2 -yv/W^inf?. (4) 

Firstly, let us review the properties of this 
system for I is fixed For this case, p and 
are canonically conjugate variable of a classical 
Hamiltonian H = ^p 2 — K y/l — p 2 cos — Flp, 
with p = = This system ex- 

hibits two qualitatively different orbits: rota- 
tion and libration . In the rotation regime, 
9 increases (or decreases) monotonically and 
p oscillates with small amplitude. In the li- 
bration regime, 9 and p oscillates around the 
equilibrium point P e {p e , 9 e ). The separatrix be- 
tween the two regimes determined by the sad- 
dle point P s (p s ,9 s ) with H a = H(p s ,0 s ). For 
large energy (H > H s ) the orbit is rotation, 
whereas for small energy (H < H s ) the orbit is 
libration. The equilibrium point and the sad- 
dle point are obtained by equating the right 
hand sides of Eqs. (||) and (||) to zero, yielding 
p e « %,9 e = and p s « ^,9 S = n. 

In the rotation regime, since 9 is increasing 
(or decreasing) monotonically with the angular 
frequency proportion to \J KC , from equation 
, we find that j is a fast oscillation current 
with frequency proportion to V KC and the up- 
per bound of the curre nt, i.e., the so-called crit- 
ical current j c w ^/l — p 2 where p a is the 
average of p. 

In the liberation regime, let p = p e + 
dp, considering C ^> K , from Eqs. (||) 
and (H]) we can obtain a pendulum equation: 
9 + KC^/T^ "pfsin(0) = 0. Then, we ob- 
tain 9 = 9 p sm.(u}t), 5p — ^r- cos(wt) in 

which uj = KC{\ — p 2 ) and P is deter- 
mined by the energy of the pendulum. So, 

one gets dp = 9 P \j %{1 - P 2 e ) ~ « 1(T 2 . 

For this case, the Josephson current is j = 



— ^\Jl—p 2 sm\9 p sin(wi)]. The critical current 

is jc = f V 1 -Pi sin(|6» p |) or j c = f^l-p 2 e 
when \0 p \ > tt/2. 

From the above discussion, for a fixed I one 
can only find the ac current in BJJ. These prop- 
erties exhibit the analog of the dc voltage case 
in SJJ. 

As it has been suggested in Ref. || , a dc cur- 
rent can be induced by moving the laser barrier 
across the trap with a constant velocity V = gi, 
and exhibit the analog of the critical behavior 
in SJJ's. Now the question arises: can this sys- 
tem exhibit some new properties which can not 
find in SJJ? To answer this question, let us re- 
view the case for SJJ. In SJJ, we know that the 
densities of cooper pair in two side of the junc- 
tions are equal to each other (if the materials 
are the same) and should hardly change when 
the dc current exists. This feature interpreting 
to the BJJ's is that the change in the relative 
population is very small. To see this feature 
clearly, let us assume the initial relative popu- 
lation p(Q) = 0, since the change in the relative 
population is very small, the Eqs. (|2|) and (||) 
can be approximated to 

p = -Ksm(0) (5) 

9 = -FVt + Cp. (6) 

Differentiating the second of these equations 
and replacing the first, we have 9 + KC sin 9 = 
—FV, which is a driven pendulum equa- 
tion. The first integral of the equation gives 
the energy of this pendulum E — \9 — 
KCcos9+FV9, and the constant E = \9{Q) 2 - 
KC cos0[(O)]+FV0(O). We define the potential 

U{0) = -KC cos 9 + FV9. (7) 

When V < 0.559^m/s, i.e., FV/{KC) < 1, 
it is a washboard potential, its local maximum 
appears at 9 — (2m — l)n + arcsin/c, and the 
local minimum appears at 9 = 2mn — arcsinfc, 
where k = and m is an integer. We re- 
strict the initial 9 in interval [— 7r,7r], then the 
motion of 9 is characterized by the local max- 
imum U c = KCy/l - k 2 + Fl/[arcsin k — 7rl. 
If E < U c , the motion of 9 is an oscillation 
in the interval [9 meix , 9 min ] with the frequency 
u ~ \/KC, where 9 max and 9 m i n is the so- 
lution of E = U{9). If E > U c , the mo- 
tion of 9 is a rotation, i,e., 9 decreases mono- 
tonically, and 9 = V / 2(E + KC cos 9 - FV9). 
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When V > 0.559/im/s, i.e., FVj(KC) > 1, the 
potential is a titled-step potential, 9 is also a 
rotation with 9 ~ —FVt. From (pi), we know it 
is an ac current when 9 is a rotation. 

To investigate the case when 9 is in oscillation 
regime, let p = p e (t) + Sp where p e {t) = ^-f-t, 

from @, we obtain Sp ps ^, so, we know that 
<5p is a small term with <5p ~ ps 0.02. Then 

from (|]) we can obtain a dc current (to zero 
order of Sp) 

FV 

* = 2C- (8) 

From the above discussion, when 9 is in oscil- 
lation regime, the current is a dc current with 
the relative population increasing (p ps p e (t)); 
when 9 is in rotation regime, the current is an 
ac one. The motion of 9 is determined by the 
initial conditions and the driven force If 

sets the initial conditions: 9(0) = 9(0) = 0, one 
can obtain a critical velocity V c = 0.406/im/s. 
For V < V c , then E < U c , it is a dc current, 
but for V > V c ,thcn E > U c , the current is in 
ac regime. This is a close analog of the critical 
behavior in SJJ. 

However, for the BJJ's, the dc current must 
lead to the change of the relative population. 
This feature will give rise to new properties. 
Taking the change of the relative population 
into account, we obtain a nonrigid driven pen- 
dulum equation 

9 + KC^\-p 2 sin 9 = -FV. (9) 

In analog to the above discussion, we get the 
potential 

U(p, 9) = -KCy/l - p 2 cos 9 + FV9. (10) 

This is also a tilted washboard potential for 
V < 0.559/im/s. The new feature is the 
local maximum U c (p) will decrease with \p\ 
increasing. If the motion of 9 is an oscillation at 
the beginning (U c (p) > E), we will firstly find 
a dc current with U c (p) decreasing. The oscil- 
lation of 9 will not keep when U c (p) < E, at 
this moment 9 will start rotating, and then the 
current will transition to an ac one. This means 
the dc current does not keep for all the time, 
i.e., the dc current has a lifetime. So one can 
have a straightforward definition of lifetime of 
a dc current as: the lifetime of the dc current t c 



is the time at which E > U c (p(r c )). This def- 
inition is consistent with the definition of the 
critical velocity in SJJ's. 

The voltage-current characteristic is the most 
important property in SJJ's. But in BJJ, the 
important physical quantity should be the rel- 
ative population p which can be directly de- 
tected. So, the critical behavior should be char- 
acterized by the change of p in BJJ. When 9 
is in oscillation regime, from the above dis- 
cussion the relative population is increasing. 
When 9 is rotating, if 9 < (arcsinfc — n), 
9 ps y/2(E + KC cos 9 - FV9) , Sp is still a 

small term with Sp max ~ \[% ~ 0.02, hence 
the relative population is still increasing and 
p ps p e (t) (in the zero order of Sp). But if 
9 > (arcsinfc — n), 9 will decrease monotoni- 
cally, 9 ps -FVt, then 6 ps -\F Vt 2 , the in- 
tegral of the current: J t °° K sin(^FVt 2 )dt ps 
0.0175 where to is a finite time, so p should 
hardly increase. This means that when the 
motion of 9 changes from oscillation to rota- 
tion, the relative population will still increase 
until 9 > (arcsinfc — tt). So, the more accurate 
definition of the lifetime of the dc current in 
BJJ should be: the lifetime t v is the time when 
@( t p) ^ (arcsin fc — 7r), within this time the rela- 
tive population is increasing, but after this time 
the relative population will hardly increase and 
keep on average fixed. 

In Fig. 1, we show the phase diagram relat- 
ing to the critical onset in the parameter space 
of t and the laser velocity V for the initial con- 
ditions: 1(0) = 0, p(0) = and 9(0) = 0. The 
solid line is the lifetime t p (V), the dashed line 
is the lifetime t c (V). The step structure implies 
that the transitions occur in different cycle of 
the oscillation. The abruptly increase is due to 
that the transition occurs near the peak of the 
potential where 9 ps 0. Fig. 2. plots the relative 
population after the dc current is destroyed, the 
crosses are the relative population p at t = 2s 
which are obtained by integrating the Eqs @j 
and (^), the solid line is the theoretical result 
given by p — p e (r p ), which shows the theoreti- 
cal estimation is consistent with the numerical 
simulation. 

From the above discussion, we know that the 
lifetime of a dc current is also determined by 
the initial conditions. One knows that the ini- 
tial relative population p(0) must be very close 
to the equilibrium p e to obtain a dc current, 
but the initial relative phase can be various. In 
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Fig. 3, we plot the lifetime of the dc current for 
different initial relative phase 0(0), where we let 
the initial population p(0) — and 1(0) = 0. 

Concerning a possible realization of the phe- 
nomenon described in this work, one should 
consider the influence on the initial acceler- 
ations. We choose the initial conditions as: 
1(0) = 0, p(0) = 0,6(0) = and V(0) = 0. 
Let the barrier starts moving with a constant 
initial acceleration. When its velocity reaches 
the value Vq , we stop accelerating and then keep 
the velocity. What we concern is the change on 
the lifetime of dc current caused by the different 
initial acceleration. Fig. 4(a) shows the lifetime 
for different initial acceleration, where the life- 
time is the time duration of a dc current after 
the acceleration. The solid line is for Vq = 0.5 
/im/s and the dashed line is for Vq — 0.4 \imj s. 
We find that for a sudden change of the initial 
barrier, the influence of the initial accelerations 
can be negligible, whereas for slow acceleration 
process the lifetime changes abruptly. The rea- 
son can be given by the following analysis: For 
a sudden acceleration, the time for accelerating 
is very short. In this time period , the change 
of 9 is small, so does the change of lifetime. On 
the contrary, in a slowly acceleration process, 
the 9 changes greatly (see Fig. 4(b)). Because 
the lifetime is very sensitive to the initial condi- 
tions, so the lifetime will change dramatically. 

We note that the conserved conditions im- 
plies that I must be less than y w 1.16/xm to 
ensure p e < 1. In another aspect, to ensure the 
number of the condensed atoms in one well be- 
yond the minimum threshold, / has to be less 
than this value. On the other hand, although 
the Eqs. (§) and (||) is obtained by solving 
variationally the GPE ([!]) under the approxi- 
mation p < < 1 , it still remains a good approx- 
imation even for p m 0.4 So, the discussion 
in this paper is accurate at least under these 
constraints. 
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FIGURES CAPTION 

Fig. 1. The phase diagram relating to the 
critical onset in the parameter space of t and 
the laser velocity V for the initial conditions: 
1(0) = 0, p(0) = and 0(0) = 0. The solid line 
is the lifetime t p (V), the dashed line is t c (V). 

Fig. 2. The relative population after the dc 
current is destroyed. 

Fig. 3. The lifetime of the dc current for dif- 
ferent initial 9. The initial relative population 
isp(0) = and Z(0) = 0. 

Fig. 4(a). The lifetime for different initial 
accelerations. The solid line is for Vq = 0.5 
/im/s and the dashed line is for Vq = 0.4 /im/s. 
(b). The relative phase 9 when the velocity 
reaches Vq. 
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